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Resolved Millimeter Emission from the HD 15115 Debris Disk 

Meredith A. MacGregor^, David J. Wilner^, Sean M. Andrews^, A. Meredith Hughes^ 

ABSTRACT 

We have used the Submillimeter Array (SMA) to make 1.3 millimeter obser¬ 
vations of the debris disk surrounding HD 15115, an F-type star with a putative 
membership in the f3 Pictoris moving group. This nearly edge-on debris disk 
shows an extreme asymmetry in optical scattered light, with an extent almost 
two times larger to the west of the star than to the east (originally dubbed the 
“Blue Needle”). The SMA observations reveal resolved emission that we model 
as a circumstellar belt of thermal dust emission. This belt extends to a radius of 
~ 110 AU, coincident with the break in the scattered light profile convincingly 
seen on the western side of the disk. This outer edge location is consistent with 
the presence of an underlying population of dust-producing planetesimals under¬ 
going a collisional cascade, as hypothesized in “birth ring” theory. In addition, 
the millimeter emission shows a ~ 3cr feature aligned with the asymmetric west¬ 
ern extension of the scattered light disk. If this millimeter extension is real, then 
mechanisms for asymmetry that affect only small grains, such as interactions with 
interstellar gas, are disfavored. This tentative feature might be explained by sec¬ 
ular perturbations to grain orbits introduced by neutral gas drag, as previously 
invoked to explain asymmetric morphologies of other, similar debris disks. 

Subject headings: circumstellar matter — stars: individual (HD 15115) — sub¬ 
millimeter: planetary systems 


1. Introduction 


Nearly a hundred dusty debris disks around nea rby stars have been spatially resolved 
at one or more wavelengths (jMatthews et ah 2014). The bulk radial structure of these 
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disks is generally well explained by the presence of a localized belt of planetesinials, or 
“birth ring,” where sma ller and smaller dust grains are produced in a catas trophic collisional 
cascade and dispersed fjStrubbe fc Chiand l2006l: lAugereau fc Beustll2006l ). However, many 


disks exhibit additional substructure such as brightness asymmetries, offsets, warps, and 
clumps that cannot be explained by the steady-state collisional models assumed in this 
framework. Numerous mechanisms have been proposed to explain morphologies that depart 
from axial syr nmetry, including i nteractions with interstellar gas, planetary r esonances, and 
stellar flybys f Pebes et ah 2009 : Maness et ah 2009 . 2008 : Kalas et ah 2007 . and references 
therein). Still, most of these structures have been detected only at optical and infrared 
wavelengths, probing a population of small grains that react strongly to stellar radiation 
and winds. Additional observational constraints, particularly at millimeter wavelengths, are 
needed to access larger grains with dynamics more similar to the parent colliding bodies, to 
solidify our understanding of the mechanisms that shape debris disk structure. 


HD 15115 is an F2V star at 45 ± 1 pc flvan Leeuwedl2007l) whose sp ace motions sug¬ 
gest membership in the young f3 Pictoris moving group (iMoor et al.l 120111) . which includes 
the well-studied (3 Pic an d AU Mic debris d isks. Although the case for membership in this 
group is not iro nclad fseelDebes et al.l 120081) . recent estimates suggest an age of 21 ± 4 Myr 


for these stars ( Binks fc Jeffries! 2014 ). An infrare d excess from HD 15115 indicating orbit¬ 


ing dust was hrst noted from IRAS observations fiSilverstonel l2000l) . Subsequent scattered 


cumstellar disk ( 

Ka 

as et al. 

2007; 

Debes et al. 

2008; 

Rodieras et al. 

2012; 

Schneider et al. 

2014: 

Mazover et al. 

2014 

). This disk shows an extreme asymmetry in optical scattered light 


(iKalas et al.l 120071) : the east side of the disk extends to ~ 7" (315 AU), while the west side 


reaches > 12" (> 550 AU). Moreo ver, there ar e indic ations that this asymmetry has roots 
at smaller scales. As described bv iDebes et al.l (l2008l) . the optical surface brightness drops 
steadily on the east side from ~ 1" (< 45 AU), while the western side appears to flatten 
interior to ~ 2" (~ 90 AU). In “birth ring” theory, this flattening corresponds to the outer 
edge of the planetesimal belt where smaller grains are created in collisions and are subse¬ 
quently launched by stellar forces into wider bound and unbound orbits. The presence of a 
ring stru cture is supported by a subtle concavity visible in scattered light to the north of 
the star fISchneider et al .1120 141) . and t his brighter side of a ring becomes clear in recent high 
contrast images (IMazover et al.ll2014l) . 


A 3a detection of 850 pm continuum emission from the HD 15115 system using the 
James Clerk Maxwell Telescope/SCUBA suggested the presence of a reservoir of large dust 
grains in t he disk, with low temperatu re (T^ = 62 K) and low mass (~ 0.047 M 0 ) indicative 
of debris ( Williams fc Andrews! I 2 OO 6 ). Subsequent mapping observations with SCUBA-2 


obtained a consistent but higher flux density (8.5 ± 1.2 mJy vs. 4.9 ± 1.6 mJy) and did 
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not resolve any strncture with a 14" (FWHM) beam fjPanie et al.ll2013f) . In light of the 
asymmetric morphology of the HD 15115 disk seen in scattered light, we have used the 
Submillimeter Array (SMA)0 to resolve the millimeter emission from the disk, aimed at 
tracing the distribution of underlying planetesimals. These observations reveal a planetesimal 
belt together with tentative evidence for an asymmetric extension of millimeter emission to 
the west, coincident with the extended feature seen in optical scattered light. We quantify 
the detected millimeter structure and discuss the plausibility of proposed mechanisms for 
the asymmetry. 


2. Observations 


We observed HD 15115 in fall 2013 with the SMA (IHo et al.l |2004J) on Manna Kea, 
Hawaii at a wavelength of 1.3 mm using both the compact and extended conhgurations 
of the array. Table [T] summarizes the essentials of these observations, including the dates, 
baseline lengths, and atmospheric opacity. With the extended baselines, these observations 
probe angular scales down to < 1". While typically only six (or fewer) of the eight array 
antennas were available during the hve tracks, the weather conditions were generally very 
good for observations at this wavelength. The total bandwidth available was 8 GHz derived 
from two sidebands spanning ±4 to 8 GHz from the local oscillator (LO) frequency. The 
phase center was located at a = 02^26“16!24, <5 = +06°17'33"19 (J2000), corresponding 
to the position of the star uncorrected for its proper motion of (86.31,-49.97) mas yr“^ 
(Ivan LeeuwenI 120071) . The ~ 54" (FWHM) held of view is set by the primary beam size of 
the 6-m diameter array antennas. 


The data from each track were calibrated independently using the IDL-based MIR 
software package. Time-dependent complex gains were determined using observations of 
two nearby quasars, J0224-I-069 (0?8 away) and J0238-I-166 (10?8 away), interleaved with 
observations of HD 15115 in a 12 minute cycle for the extended tracks and a 16 minute cycle 
for the compact tracks. The passband shape was calibrated using available bright sources, 
mainly 3G84 or 3G454.3. Observations of Uranus during each track were used to derive the 
absolute hux scale with an estimated accuracy of ~ 10%. Imaging and deconvolution were 
performed with standard routines in the MIRIAD software package. A variety of visibility 
weighting schemes were used to explore compromises in imaging between higher angular 
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Academica Sinica Institue of Astronomy and Astrophysics and is funded by the Smithsonian Institution and 
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resolution and better surface brightness sensitivity. 


3. Results and Analysis 


3.1. 1.3 mm Dust Continuum Emission 


Figure 1 shows a contour image of the 1.3 millimeter emission overlaid on a Hubble 
Space Te/escojpe/Ad vanced Came r a for S urveys coronographic image of optical scattered light 
(F606W hlter) from iKalas et al.l (120071) . The synthesized beam size for this 1.3 millimeter 
image, obtained with natural weighting, is 3'.'0 x 2'.'1 (135 x 95 AU), position angle 73°. The 
rms noise is 0.30 mJy beam“^, and the peak signal-to-noise ratio is about 5. The red star 
symbol indicates the stellar position, corrected for proper motion, i.e. offset from the phase 
center by (1'.'3, —0/7). This position coincides very closely with the emission peak, well 
within the combined uncertainties of the absolute astrometry and modest signal-to-noise. 
The image reveals a narrow band of 1.3 millimeter continuum emission that extends from 
the east-southeast to west-northwest, aligned closely with the optical scattered light disk 
orientation. The western extent of the millimeter emission appears to be greater than the 
eastern extent, although the signihcance of this elongation is not high. 



Fig. 1.— SMA image of the 1.3 millimeter continuum emission from HD 15115, overlaid on an 
image of optical scattered light from the Hubble Space Telescope (jKalas et al.l 120071 ). The contour 
levels are in steps of 2 x 0.3 mJy, the rms noise level. The ellipse in the lower left corner indicates 
the 3/0 X 2'.'! (FWHM) synthesized beam size. The star symbol marks the position of the stellar 
photosphere, corrected for proper motion. 
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3.2. Disk Modeling 


In order to characterize th e millim eter emission from HD 15115, we used the procedure 
described by MacGregor et al.l (120131) that employs Markov Chain Monte Carlo (MCMC) 
methods to £t a simple parametric model to the observed visibilities. We assume that 
the emission arises from a geometrically thin, axisymmetric belt characterized by a surface 
brightness prohle of the form Iu{r) oc for < r < i?out- The belt emission normalization 
is dehned by a total flux density Fbeit = / Rdfl, and the belt center is given by offsets relative 
to the pointing center {Aa, A(f}. The various scattered light observations show that the disk 
is nearly edge-on to the line of sight. For our models, we adopt an incli nation angle of 87 ° 
and an orientation on the sky described by a position angle of 278?5±0.5 (IKalas et al.ll2007l) . 
Small variations in these values have no material impact on the results. This initial model 
does not attempt to account for any asymmetric structure in the emission. To address the 
western extension, we also made models that included an unresolved source with flux density 
Fj-es and position Ar, dehned as an offset from the belt center along the position angle of 
the disk major axis. 

For a given parameter set, we compute synthetic visibilities corresponding to the SMA 
observations and compare directly to the data with a value (the sum of real and imaginary 
components over all spatial frequencies). By htting the visibility data directly, we are not 
sensitive to the non-linear effects of deconvolution, and take full advantage of the complete 
range of spatial frequencies sampled by the observations. We assumed uniform priors for 
all parameters, with reasonable bounds imposed to ensure that the model was well-dehned: 
Tbeit A 0 and 0 < i?in < Rout- In addition, we constrained the two parameters {Aa, A(5} 
that describe the belt center to be within 0'.'5 of the offsets predicted from the stellar proper 
motion; this generously accomodates the uncertainties in the proper motion and the absolute 
astrometry of the observations. The ht quality is characterized by a likelihood metric, £, 
determined from the values (l n£ = — y^/ 2). We make u se of the affine-invariant ensem¬ 
ble MCM C sampler proposed bvICoodrnan fc Wearel (120lOh and implemented effectively in 
Python bv iForeman-Mackev et ahl (l2013l) . With this MCMC approach, we can characterize 
efficiently the multidimensional parameter space of this simple model and determine pos¬ 
terior probability distribution functions for each parameter by marginalizing over all other 
parameters in turn. 


3.3. Results of Model Fits 

Figure 12] shows the output of ~ 10'^ MCMC trials and Table [2] lists the best-£t model 
parameter values and their 68% uncertainties determined from the marginalized posterior 
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probability distributions for the full model (axisymmetric belt plus western extension). Fig¬ 
ure 3 shows comparisons between the data and the best-£t models in the image plane, 
including the imaged residuals in the rightmost panels. The simple symmetric disk model 
in the upper panels of Figure 3 reproduces the bulk of the millimeter emission well, but the 
residual image shows a ~ 3(T feature coincident with the western extension of the scattered 
light along the disk axis, and this model yielded a reduced = 3.67 (153,032 independent 
data points, 6 free parameters). The departure from unity reflects the poor ht of the model, 
in part due to the presence of the residual western emission. By including additional pa¬ 
rameters in the model that account for this residual feature (F).es and Ar), the model fit was 
more satisfactory, as can be seen in the lower panels of Figure 3, with reduced = 1.88. 


The total flux density of the best £t model is constrained to be Fkeit = 2.6log niJy. If 
we extrapolate the available submillimeter measurements in th e literature using th e typical 
spectral index of ~ 2.65 for debris disks at these wavelengths flGaspar et ahl 120121) . we hnd 
cl ose agreement with the SMA value. In particular, extrapolating the SCUBA-2 observation 
ofiPanic et ahl (12013h to 1.3 mm gives 2.8 ± 0.4 mJy, while the older SCUBA observation of 


Williams fc Andrews! (120061) gives 1.6 ± 0.5 mJy. Civen this consistency, it seems likely that 


the SMA observations detect the full disk emission and do not miss any signihcant, more 
extended, millimeter flux. 


The location of the outer edge of the modelled millimeter emission belt, i?out = IIOI 22 
can be compared to expectations based on the scattered light ra dial surface b r ightne ss pro- 
hle and “birth ring” theory. The E.ST imaging of HD 15115 bv iDebes et ahl (120081) shows 
that the scattered light prohle steepens on the western side of the dis k, falling as 7-- 3-56ifco.o6 


beyond 1'.'8 (81 AU). The ground based scattered light observations of iRodigas et ahl (120121) 
show a similar break in the western surface brightness prohle near ~ 1.8 — 2" (~ 80 — 90 
AU). The model ht to the outer edge of the millimeter belt is consistent with both of these 
determinations, within the substantial uncertainties. However, in the context of theory, it 
remains unclear why the eastern side of the scattered light disk shows no clear break in 
surface brightness slope. 


The observations of iRodigas et ahl (120121) additionally show a drop in the scattered 
light surface brightness of the disk’s western side interior to 1" (45 AU), presumably m ark¬ 
ing the inner edge o f the belt. The higher contrast images of ISchneider et al.1 (120141) and 
Mazover et ahl (12014J 1 conhrm this feature and suggest at least a partial clearing with i n this 
radius. Analysis of the spectral energy distribution of the system by iMoor et ahl (120111) 
yields a best-ht model with an inner warm component and an outer cold ring at a minimum 
radius of 42 ± 2 AU. Both of these determinations of the inner disk edge match the (poorly 
constrained) best-ht inner radius of the millimeter belt, Rin = 43 ± 28 AU. 
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Despite the modest signal-to-noise ratio of the millimeter observations, this model htting 
exercise yields constraints on the inner and outer radii of the millimeter emission belt that 
are consistent with previous observations and compatible with the “birth ring” model of 
debris disks. 


Not surprisingly, given the limits of the resolution and sensitivity of the SMA data, 
model htting does not provide a strong constraint on the power-law gradient of the millimeter 
radial emission prohle, x = —O.TS^q^. The well known degene racy between the emission 
gradient, x, and disk outer radius, i?out, (e.g. iMundv et ahlll9961 ) is apparent in Figure [2] as 
a slope introduced in the lower contours resulting from that pair of model parameters. For 
very negative values of x, the contours spread to span a wide range of potential outer radii, 
as the width and extent of the emission belt become difficult to constrain due to the low 
brightness in the outer regions. Similarly, the inner radius of the disk is poorly constrained 
for positive gradients. If we assume that the emitting dust is in radiative equilibrium with 
stellar heating, leading to a temperature gradient described by T oc then taking the 

best ht exponent at face value implies a relatively shallow surface density prohle, S oc 
(albeit with large uncertainty). 


4. Discussion 


Millimete r imaging of HD 15 115, like its su spected sister stars in the 0 Pictoris movin g 
group, (3 Pic ( Wilner et ah 2011 1 and AU Mic f Wilner et ah 2012 : MacGregor et ah 2013), 
shows a resolved belt of emission with outer edge coincident with a previously observed break 
in scattered light, consistent with an underlying “birth ring” of colliding planetesimals. For 
all of these stars, resolved multi-wavelength datasets suggest planetesimal collisions within a 
belt produce grains with a spectrum of sizes, and the ehects of size-dependent dust dynamics 
produce the compact millimeter emission and an extended scattered light halo. HD 181327, 
a fourth member of the (3 Pictoris moving group, a lso shows a belt of c old dust seen in 
Herschel PACS 70 pm and ATCA 3.2 mm images (jPebreton et ah 2012 1. though higher 
resolution millimeter images are still needed to determine its underlying structure. 

The detection of a ~ 3cr residual along the western extent of the HD 15115 disk, aligned 
with the asymmetry seen in scattered light from small grains, provides tantalizing evidence 
that the distribution of larger grains in this system may be asymmetric as well. Similar 


asymmetric debris disk morphologies in scattered light have been obser ved in other 


inclined systems where the geometry favo rs detection, nota bly HD 32297 fICurrie et al 


Debes et al. 

2009 

2009 

), as well as 


Boccaletti et ah 


2012) and HD 61005 ( Hines et al.l 120071: iManess et al. 


lighly 


2012 


2 OO 9 I) . as well as (3 Pic ( Kalas fc Jewitt 1995 ). While all of these disks show unique mor- 


















































phological details, similar physical mechanisms may be responsible for the overall shaping 
of the asymmetric appearance of the debris. However, it is unclear if the asymmetries can 
be clearly traced back to the dust-producing parent bodies. The constraints from millimeter 
observations so far fail to provide a co nsistent picture. H D 32297 was observed at 1.3 mil¬ 
limeters in the early days of CARMA flManess et al.l 120081) . and the millimeter image shows 
a hint of asymmetry, albe it at low signal-to-n oise. On the other hand, SMA 1.3 millimeter 
observations of HD 61005 flRicarte et al.ll2013h provide no evidence that the millimeter grains 
are present in the spectacular swept-back wings seen in scattered light. 

The most common mechanisms invoked to explain the scattered light asymmetries in 
these highly inclined debris disks involve interactions with the local interstellar medium 
(ISM). Ram pressure from interstellar gas can remove bound and unb ound grains from a 


disk directly, or by inducing outflows of disk gas that entrain the grains flManess et al.ll2009 


Debes et al.l 120091) . Additionally, neutral gas can introduce secular perturbations to the or- 
bits of bound yains , producing signihcant effects on disk morphology on orbital timescales 
flManess et al.ll2009l) . Alternatively planets within the system can hav e dust-trappi ng reso¬ 
nances that create clumpy, wavelength-dependent grain distributions flWvattI 120061) . Look¬ 
ing outside of the system entirely, ext ernal perturbation by a close stellar flyby can also 
shape the distribution of disk grains flKalas et al.l 120071) . An increase in the dust grain 
scattering cross section or a l ocal density enhancement could produce a local increase of 
collisions flMazover et al.ll2014l) . Finally, azimuthal asymmetries could be generated through 
the “ photoelectric instability ” effect, if the disk has a high enough total dust to gas mass 
ratio flLvra fc Knchnerll2013l) . Each of these mechanisms affects different sized grains pref¬ 
erentially and works under specihc physical conditions. In the following, we examine the 
plausibility of each of these mechanisms for creating asymmetric structures in the HD 15115 
disk and similar systems. 


4.1. Interactions with the ISM 


Because the p roper motion of HD 15115 is nearly parallel to the disk major axis, 
Debes et al.l fl2009l) suggest that the dominant disk asymmetry in scattered light arises from 
an ongoing interaction with the local ISM. In this picture, the eastern side of the disk plows 
head hrst into the ISM, causing that side to become truncated, while small grains are blown 
out to the west. Motiva t ed by the swept back morphologies of HD 32297, HD 15115, and 
HD 61005, iDebes et ahl fl2009l) construct a model for a disk interacting with the gaseous 
ISM, based on the gas drag felt by meteorites entering the Earth’s atmosphere. This results 
in a scaling law for the radial distance (Rdeflect) at which unbound grains are signihcantly 
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perturbed from their original circumstellar orbits: 


R 


deflect 


= 286 


1.67 X 10 g cm ^ 

Pgas 


20 km s 

'^rel 


-1 




( 1 ) 


where pgas is the ISM gas density and Urei is the relative cloud-disk velocity. By applying the 
dehnition of f3, the ratio of the force of radiation pressure to the gravitational force exerted by 
the star, Equation (1) can be recast into a scaling relation for the ISM gas density (ugas) and 
disk relative velocity (urei) needed t o perturb un^bound grains from their orbits. Assuming 
perfectly absorbing spherical grains (IBurns et al.lll979[) 


(3 = 0.574 


Lq 


Me 


-1 


1 g cm ^ 

Pdust 


1 fim 
?^dust 


( 2 ) 


which yields the relation 


n. 


gas 


200 cm 3 / V 30 km s 






> 

rsj 


0.6262 


Af* \ f Odust 


0.95 Mq J \0.1 pm 


X 


Pdust 


2 g cm 



70 AU 


( 3 ) 


Interpreting HST images of HD 61005. iManess et al.l (1200911 discuss a second model for 
ISM interaction in which an interstellar gas cloud removes both bound and unbound grains 
from a disk due to ram pressure stripping. In this scenario, the drag force on a grain must 
be compara ble to or exceed the gravitational force binding the grain to the star. For the 
bound case, iManess et al.l (120091) show that this leads to the following scaling relation for 
the cloud density (ugas) and the relative cloud-disk velocity (again Urei) 


/ ^gas \ f ^rel \ ^ ( Af* \ / Odust \ 

V200 cm-sy V30 km s"V ~ VO-95 Me) VO-l pm/ 

X 

\2 g cm V ^ 


( 4 ) 


A similar relation is obtained for the unbound case with an additional constant factor ac¬ 
counting for the relative distance traveled by a grain parallel and perpendicular to the disk 
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midplane. Although formulated somewhat differently from iDebes et al.l (120091) . the result¬ 
ing expression is functionally identical with the exception of a constant factor (compare 
equations [3] and 0]) . 

While ram pressure stripping of disk grains by the IS M seems like an attractive ex pla- 
nation for the asymmetry seen in systems like HD 15115, iMarzari fc Thebaultl (1201111 use 
numerical modeling to show that even grains just above the disk blow-out size (1 — 10 /rm) 
are minimally affected by ISM interactions and are lost from the disk before they acquire 
any significant inclination. The millimeter grains, with much smaller values of /3, should feel 
little to no effect from ISM gas. Thus, this argument has difficu lty explaining an asymmetry 
in the millimeter disk emission. Moreover, iManess et al.l (120091) point out that the ISM den¬ 
sities required are characteristic of cold clouds (u ~ 50 cm“^, T ~ 20 K) that occupy only a 
very small volumetric filling factor within the Local Bubble. 


Maness et al.l (I2009|l also propose that disk gas could undergo ram pressure stripping 


by the ISM, and disk grains are consequently swept away when they become entrained in 
this outflowing gas. In this scenario, sufficient disk gas is required and the grains must 
be entrained effectively, neith er of which is a ssured. However, measurements of the NA I 
doublet towards HD 15115 bv iRedfieldl (120071) indicate that the column density for this disk 
is ~ 5 times greater than towards (3 Pic; he estimates that the upper limit on the total gas 
mass in the circumstellar disk is O.3M0. 

In order for direct ram pressure stripping of disk grains to occur, t he force of the ISM on 
the grains must be comparable to the stellar gravitational force. Yet, iManess et al.l (1200911 
propose that in the case where this condition is not met, neutral gas can still affect disk 
morphology over timescales of ~ 10^ — 10“^ years by introducing secular perturbations to the 
orbits of bound grains. A similar mechanism has been proposed as the removal mechanism 
for dust from our Solar System at 20 — 100 AU (jSchererl 120001) . Like drag from the Solar 
wind, neutral gas drag involves the transfer of momentum from incident gas particles to 
grain surfaces. However, unlike the Solar wind, interstellar gas drag tends to increase grain 
eccentricities and semi-major axes, and can eventually unbind grains from the system. The 
normalcy of the interstellar densities, velocities, and cloud sizes required by this neutral gas 
drag model make it a promising mechanism for producing the structures in this family of 
asymmetrical debris disks. Furthermore, the presence of detectable gas in the disk lends 
additional credibility to the argument that neutral gas drag may play a role in shaping the 
disk morphology. 
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4.2. Interactions with Other Perturbers 


An alternative mechanism to ISM interactions that can produce disk asymmetries in¬ 
volves planet induced resonances. Parent planetesimals locked in resonance with orbiting 
planets can produce large grains that stay in resonance, or, if collisions are unirn p ortant. 


grains can drift into resonances due to Poynting-Robertson drag flKrivov et ahll2007() . IWvatt 


(1200611 predict that in the former case, an asymmetry can be present at short and long wave¬ 
lengths, but absent at intermediate wavelengths. In this model, the large grains that domi¬ 
nate millimeter emission trace the parent planetesimals and thus exhibit the same clumpy, 
resonant distribution. The small grains dominant at short wavelengths are preferentially cre¬ 
ated in the high-density resonant clumps and expelled from the system on short timescales. 
The mid-sized grains traced by intermediate wavelengths remain bound, but fall out of 
resonance due to ra diation pressure and are scattered into an axisymmetric distribution. 
Maness et al.l (1200811 invoke this prediction to explain the multi-wavelength observations of 
HD 32297, which appears asymmetric at both optical and millimeter wavelengths, but sym¬ 
metric in the mid-infrared. A similar picture may be hinted at in the case of HD 15115. 
The scattered light observations are clearly asymmetric, and the SMA observations suggest 
a potential m illimeter asymmetry . However, observations at 3 — 5 pm show a symmetric disk 
morphology (IRodigas et al.l 1201211 . Detailed dynamical models of these systems are needed 
to draw hrm conclusions about this speculation. 

Global disk asymmetries could also be produced by a single or periodic stellar flyby, 
as proposed to explain the large scale s cattered ligh t asym metry of the (3 Pic debris disk 
( Kalas fc Jewitt 1995 : Kalas et ah 2001 1. Kalas et ah ( 2007 1 poi nt out that the /3 P ic moving 
group member, HIP 12545, nearby in the sky to HD 15115 f3?9 lMo6r et al.ll2006[l . could be 
involved in a past interaction, as these two stars show identical proper motion, galactic space 
motion, and heliocentric distance within the measured uncertainties. Morever, the present 
location of HIP 12545 is in the direction of the truncated eastern side of the HD 15115 


disk, as in the simulations of (iLarwood fc KalasI 1200111 flyby interactions. However, there 
are no indications other than this circumstantial evidence for any interaction. Furthermore, 
the membe r ship of HD 15115 in the [3 Pic moving group has been called into question. 
Moor et al.l (120061) us ed kinemati c argu ments to propose this group membership. The recent 
Bayesian analysis by iMalo et al.l (120131) suggests a higher probability of membership in the 
Columba Association. A better understanding of the group membership and age of HD 
15115 certainly would be helpful in unraveling the origins of the debris disk morphology. 
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4.3. Interactions within Disk Material 


_ ^Gemini Near-Infrared Coronographic Imager (NICI) data in the H and Ks bands flMazover et al, 

2 OI 4 J) resolved the HD 15115 disk and detected ansae on both sides at a radius of 1.99" (~ 90 


AU). These new observations indicate a ring-like shape with an inner cavity that appears 
symmetric, in contrast with the east-west brightness asymmetry seen in other images. If 
the inner ring is truly symmetric, an alternative o r add itional mechanism for the observed 
brightness asymmetry is possible. iMazover et al.l (120141) suggest that both a local increase 
in the dust grain scattering cross section and/or a local density enhancement of small grains 
could lead toward an increase in collisions at the location of the western brightness peak. A 
similar explanation was recently invoked t o explain the mil limeter emission clump seen in 
ALMA observations of the /? Pictoris disk ( Dent et al.l 12014 ). 


We can crudely estimate the mass of small grains needed to account entirely for the 
tentative excess millimeter emission on the western side of the disk. Optically thin dust 
emission at a temperature, Tjust, has flux F^ust ~ ^i/-B!/(^dust)Afdust/d^, where Hi, is the 
dust opacity, Hj, is the Planck function, Mdust is the dust mass, and d = 45 pc. We estimate 
a characteristic temperature for dust in this excess feature of ~ 30 K, using its approximate 
radial l ocation (^160 AU) and a ssuming radiative equilibrium of starlight with blackbody 
grains (IBackman fc Parescelll993[) . To estimate an appropriate dust opacity, we assume that 
the grain size distribution is described by a power-law between 0.1 pm and some maximum 
size Umax, n{a) oc a‘^. For grains with composition volume fractions of 0.07, 0.21, 0.42, 


and 0.30 for silicate, carbon, ice, and vacuum, respectively, iRicci et al.l (120101) calculate 
dust opacities as a function of Umax- (Note, the results are not very sensitive to other 
reasonable assumptions about composition.) Adopting a typical power-law exponent q = 
—3.5, and maximum grain size Umax = 0.1 millimeter, models yield a dust opacity, ~ 3 
cm^ g-i. Given these assumptions, the mass of small grains required to account for the 
excess millimeter flux is ~ 1 x 10^® g, or ~ 1.4 times the lunar mass. As this large mass 
indicates, a very signihcant collision would be required for small grains to account entirely for 
the western millimeter emission feature. Furthermore, such a collisional feature is transient 
and can persist for only a small fraction of the system age (I Jackson et al.ll2014l) . 


Lyra fc Kuchnerl (120131) discuss an additional mechanism by which disk gas can generate 
azimuthal asymmetries through the “photoelectric instability” effect. In this scenario, a 
dense region of dust heats the gas through photoelectric heating and creates a local pressure 
maximum, which in turn attracts additional dust. The resulting ins tability could lead to the 
development of rings, spirals, or other structures within the disk. iLvra fc Kuchnerl (120131) 
suggest that this effect will occur in disks with a total disk dust to gas mass ratio, e ~ 1. 
In the case of HD 15115, estimates of the total dust and gas mass in the disk suggest that 
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e > 0.16, in the regime where this mechanism could work. However, without further analysis, 
it is unclear whether the level of azimuthal asymmetries produced would be large enough to 
account for these new observations. 

In short, there is still no clear consensus as to which of these mechanisms is responsible 
for shaping the debris around HD 15115 and other, similar systems (notably HD 61005 and 
HD 32297). Each mechanism affects different sized grains preferentially and operates under 
specihc physical conditions. Improved constraints on the detailed morphologies of these 
systems are required to reach dehnitive conclusions. 


5. Conclusions 


We present SMA 1.3 millimeter observations that resolve the emission from the HD 
15115 debris disk, the hrst at such a long wavelength. The bulk of the millimeter emission is 
described well by a symmetric belt, consistent with the presence of a “birth ring” of collisional 
planetesimals indicated by previous observations of scattered light. The SMA observations 
also show a ~ 3(7 residual feature coincident with the western extension of the scattered light 
along the disk major axis. If real, this additional feature hints that the distribution of larger 
grains may be asymmetric as well. 


There is no consensus concerning the physical mechanism or mechanisms responsible 
for the scattered light asymmetries in HD 15115 and other, similar debris disks. Ram pres¬ 
sure from surrounding interstellar gas can strip both bound and unbound grains, or induce 
outflows of disk gas that entrain grains. However, the tentative millimeter emission asymme¬ 
try disfavors any process that affects only small grains. An alternative mechanism invokes 
neutral gas in the disk that introduces secular perturbat ions to grain orbits, producing sig- 
nihcant effects on disk morphology on orbital timescales flManess et al.ll2009l) . The presence 
of atomic gas towards HD 15115 suggests that this idea should be explored more thoroughly. 
Interactions with perturbers are another possibility. Dust-trapping resonances from orbiting 
planets may be able to explain the observed wavelength-dependent structure. Future obser¬ 
vations of HD 15115 at higher angular resolution and sensitivity with the Atacama Large 
Millimeter/Submillimeter Array are needed to resolve the detailed morphology of the mil¬ 
limeter emission and to address the processes responsible for shaping this remarkable debris 
disk. 
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Table 1. Submillimeter Array Observations of HD 15115 


Observation 

Date 

Array 

Config. 

Number of 

Antennas 

Baseline 
Lengths (m) 

LO Ereq. 
(GHz) 

HA 

Range 

225 GHz atm. 
Opacity*^ 

2013 Sep 15 

Extended 

6 

25 - 

- 174 

225.5 

-4.2,3.3 

0.08 

2013 Oct 19 

Extended 

5 

20 - 

- 182 

235.5 

-3.8,4.9 

0.28 

2013 Nov 18 

Extended 

6 

23 - 

- 174 

225.5 

-2.3,3.1 

0.10 

2013 Dec 10 

Compact 

6 

8 - 

- 53 

225.5 

-1.3,3.7 

0.20 

2013 Dec 12 

Compact 

6 

6 - 

- 53 

224.9 

-4.5,4.8 

0.11 


Note. — characteristic value for the track measured at the nearby Caltech Submillimeter 
Observatory 


Table 2. Model Parameters 


Parameter 

Description 

Best-Pit 

68% Confidence Interval 

Tbelt 

Belt flux density (mJy) 

2.56 

+0.50,-0.83 

Rin 

Belt inner radius (AU) 

43.4 

+28.3,-28.3 

Rout 

Belt outer radius (AU) 

113. 

+30.6,-21.8 

X 

Belt radial power law index 

-0.75 

+1.37,-0.87 

res 

Point source flux (mJy) 

0.84 

+0.16,-0.12 

Ar 

Point source offset (") 

-3.55 

+0.45,-0.12 

Aa 

R.A. offset of belt center (") 

1.26 

+0.07, -0.05 

A6 

Deck offset of belt center (") 

-0.78 

+0.09,-0.05 
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R,„ [AU] 



Fig. 2.— A sample of the output from a run of 10^ MCMC trials for the 6 best-fit model belt 
parameters {Rm, Rout, Fheit, x, F^es, and Ar). The diagonal panels show the ID histogram for each 
parameter marginalized over all other parameters considered in the model. For each parameter, 
the peak of each histogram is taken to be the best-fit value. The remaining panels show contour 
plots of the Icr (red) and 2cr (gray) regions for each pair of parameters, with the blue crosses 
marking the best-fit values. The inset panel in the upper right shows the resulting best-fit model 
at full resolution, pixel scale ~ (^!02 (0.9 AU). The la and 2cj regions are determined by assuming 
normally distributed errors, where the probability that a measurement has a distance less than a 
from the mean value is given by erf 
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Fig. 3.— (upper): (left) The 1.3 millimeter continuum emission from HD 15115 observed with the 
SMA, as in FigurelH (center) image of the best-fit azimuthally symmetric disk model (see Section[3] 
for a description of the modeling formalism and results), and (right) the imaged residuals from the 
symmetric model, showing a Scr feature on the western side of the disk, (lower): (left) the SMA 
image, (center) image of the best-fit asymmetric model, and (right) the imaged residuals from the 
asymmetric model, which do not show any remaining significant features. The contour levels are 
at 2a (0.6 mJy beam“^) intervals in all panels. 












